Thermodynamic Study of Excitations in a 3D Spin Liquid 



Y.K. Tsui, J. Snyder, and P. Schiffer 

Department of Physics, Pennsylvania State University, University Park, PA 16802. 



Abstract 



In order to characterize thermal excitations in a frustrated spin liquid, we 
have examined the magnetothermodynamics of a model geometrically frus- 
trated magnet. Our data demonstrate a crossover in the nature of the spin 
excitations between the spin liquid phase and the high-temperature param- 
agnetic state. The temperature dependence of both the specific heat and 
magnetization in the spin liquid phase can be fit within a simple model which 
assumes that the spin excitations have a gapped quadratic dispersion relation. 

PACS numbers: 75.30.Kz, 75.50.Ee 
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In recent years, there has been increasing interest in the properties of spin liquid states 
in which the spins are highly correlated but have strong fluctuations in the low temperature 
limit. These states are an important component of the physics in a variety of systems includ- 
ing low dimensional magnets 0, cuprate superconductors 0, and geometrically frustrated 
antiferromagnetic materials in which the topology of the lattice leads to frustration of the 
antiferromagnetic exchange interactions 0. In these geometrically frustrated magnets, the 
spins do not order except at temperatures well below the Curie- Weiss temperature (0 C w), 
and in the low temperature limit they can exhibit a variety of novel correlated spin states 
including the spin liquids in which the spins fluctuate within low-moment locally correlated 
spin clusters [|]]. These frustrated spin liquids have been the subject of extensive theoret- 
ical work H, and their existence has been confirmed experimentally in several materials 
of different frustrated topologies P~P^I- The basic thermodynamic properties of frustrated 
spin liquids are, however, only beginning to be explored |5| JT0|Jl^Jl3| , and there is neither 
a clear understanding of the nature of the thermal excitations at low temperatures nor a 
clear distinction between the expected excitations in spin liquids and those in paramagnetic 
states at higher temperatures. 

In this paper we examine the excitations in a three-dimensional frustrated spin liquid 
through measurements of the specific heat (C) and static magnetization (M) in a model 
geometrically frustrated magnet. The field dependence of C(T) suggests a crossover in the 
nature of spin excitations between the spin liquid phase and the high-temperature param- 
agnetic phase. Moreover, both C(T) and M(T) at low temperatures can be quantitatively 
described by a spin-wave-like model of excitations with a gapped quadratic dispersion rela- 
tion. The gap grows linearly with magnetic field, but appears to have a finite value even in 
the extrapolation to zero field. 

We studied the model geometrically frustrated magnet Gd 3 Ga 5 0i2 (gadolinium gallium 
garnet or GGG) 0,f|,[ll],|l^,[l^,[il| in which the magnetic Gd ions are located on two inter- 
penetrating corner-sharing triangular sublattices within the garnet structure. This geometry 



frustrates the antiferromagnetic nearest neighbor exchange interactions (|O cw | ~ 2K) [|T6|JT7 
and yields a highly unusual low-temperature phase diagram as shown in the inset to figure 
|TJ 1000. At low fields (H < 0.08 T) previous studies show that the ground state of GGG 
is a spin glass with T g ~ 0.13 K @. At fields between 0.7 and 1.4 T the ground state is a 
field-induced long-range-ordered antiferromagnetic state with Tjv(max) ~ 0.35 K |[7|,p|,p 



The ground state at intermediate fields, 0.08 ^ H ^ 0.7 T, shows no evidence of any sort 
of spin ordering upon cooling from the high-temperature paramagnetic state down to the 
lowest temperatures which have been studied (T < |9 cw |/40), and is apparently a homoge- 
neous three-dimensional spin liquid phase in which the spins are strongly fluctuating even 
at the lowest temperatures |7| JTT| . p~| . 

We measured M and C on samples of GGG cut from the same single crystal grown 
by the Czochralski method. Magnetization measurements were performed by mounting 
a GGG sample on a flexible silicon paddle and applying a magnetic field gradient (< 1 
T/m) in addition to a spatially homogeneous magnetic field parallel to the [100] direction. 
The resulting force on the magnetized sample was proportional to the magnetization, F = 
MdH/dz, and was measured from the deflection of the paddle by a capacitive technique 
|7|jl9|. The sample, the magnetometer, and the sample thermometers were immersed with 
liquid helium and thermally linked to the mixing chamber of a dilution refrigerator. As a 
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result, we could ensure the sample was in thermal equilibrium with the refrigerator during 
the measurements. Details of the experimental setup and techniques for heat capacity 
measurements can be found in Ref. [|Il[ . 

In figures [I] and |^ we plot M(H) and C(H), and in both data sets there are features 
for the boundaries of the antiferromagnetic phase (clearly visible as sharp peaks in DM/ dH 
in agreement with previous a.c. susceptibility data |7j and recent Monte Carlo results ||). 
Both M and C change monotonically with field in the spin liquid phase (0.08 ^ H < 0.7 T), 
with M(H) increasing almost linearly with field at low temperatures. We note that C(H) 
is only weakly field-dependent in the spin liquid phase at low temperatures, changing in 
magnitude by only ~ 30 % when the applied field is doubled from 0.3 to 0.6 T at T = 
0.065 K. Taking the spin 7/2 of the Gd 3+ ions to set the energy scale, the AH of 0.3 T is 
quite large, equivalent to AT ~ 0.3/iGd/fcB ~ 1-4 K, which is comparable to the exchange 
energy between spins (|O cw | ~ 2 K). This weak field dependence of C(H) is similar to that 
of another geometrically frustrated magnet with a spin liquid state, SrCr 9p Gai2_9pOi9 pif , 
and is theoretically expected for a spin liquid in which the excitations are among correlated 
spin-cluster singlets [§,|13|]. 

An unresolved issue in the study of geometrically frustrated spin liquids is whether the 
spin liquid phases can be differentiated from the high-temperature paramagnetic states since 
there are no sharp features in the temperature dependence of physical properties to indicate 
a transition to a different phase. This is the case for the spin liquid phase of GGG, but the 
specific heat data do show a distinct crossover in behavior at T cross ~ 0.235 K, where for 
all fields less than 0.7 T, the C{T) curves cross (figure 0). This crossing indicates that all 
field derivatives of the specific heat vanish at T cross , which in turn implies that dC(H)/dH 
changes sign at T cross . The field independence of the heat capacity at T cross is also evident 
in figure |2| in that the 240 mK isotherm is constant for magnetic fields within the spin liquid 
phase. The intrinsic nature of this crossover is confirmed by an inflection point in our M(T) 
data (figure |j) which is thermodynamically equivalent to a zero in dC (H) / dH (see equation 
H below). 

This crossover temperature suggests that we might draw a distinction between the nature 
of spin excitations in the spin liquid state and those in the paramagnetic state at higher 
temperatures. The change in the sign of dC (H) / dH implies that the excitations are more 
ferromagnetic-like in the low temperature limit in that they are suppressed by a field, and 
that the excitations are more antiferromagnetic-like at higher temperatures in that they are 
enhanced by a field. In other words, at high temperatures the thermal fluctuations primarily 
dissociate the antiferromagnetic correlations between the spins, while at low temperatures, 
the antiferromagnetic correlations are unaffected by thermal fluctuations so the fluctuations 
primarily reduce the alignment of the spins with the magnetic field. This crossover may 
actually be a generic feature of spin liquid states, since the same feature was recently ob- 
served by Ramirez et al. in SrCr 9p Gai2-9pOi9 These two cases suggest that a generic 
distinction can be drawn between the low temperature spin liquid and the high temperature 
paramagnetic states in these systems, at least in the nature of the spin excitations in the 
presence of a magnetic field. 

To further investigate the nature of the low temperature excitations in the spin liquid 
phase, we closely examined our C(T) data for T< T cross within the spin liquid regime. Since 
C(H) shows ferromagnetic-like behavior (i.e. dC(H)/dH < 0), we fit the low temperature 
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C(T) data to a spin- wave-like model with the dispersion relation, Hu = A + Dk 2 , where D is 
the spin stiffness constant and A is an energy gap. Note that while this dispersion relation 
has the same form as that for ferromagnetic spin waves, in fact we are only assuming that 
the excitations are non-interacting bosons with a gapped quadratic dispersion relation. In 
this model, we calculate the specific heat as 

_ 4 /2t3/2 r gVj - x exp{x) , m 
model1 J ~47rW 2 4 (exp(x)-l) 2 ' [) 



where x = A/k^T. As shown in figure [|, equation [I] fits the measured C(T) data well [p22 
and we obtain D ~ 2 x 1CT 43 Jm 2 which is almost independent of field and A which is linear 
in field as shown in the inset to figure H. Moreover, since 




'dhf 



(2) 



H 



we can get the temperature dependence of M by integrating the expression for the specific 
heat C, 

M(T)~ £ dT> £' ^p-dT" . (3) 

We perform this integral numerically, substituting equation |T] in equation |3|, and we find 
that the numerical form for the magnetization can be described by 

M(T) ~ M (l — AT 3 ), (4) 

in the temperature range where we fit C(T) (60 ^ T < 200 mK). As shown in figure 
|5], equation [| does fit the measured M(T) very well at low temperatures, lending further 
credence to the model. 

Since the model apparently describes the measured values of C(T) and M(T) quite well, 
we can examine the obtained fit parameters to obtain some understanding of the nature 
of the excitations. The fitted values of the spin stiffness parameter, D, vary by less than 
10% in our range of fields. If we take the usual assumption for a ferromagnet that the 
ferromagnetic exchange energy J = D/2Sa 2 where S = 7/2 for the Gd ions and a = 3 A is 
the nearest neighbor Gd separation, we obtain J ~ 0.02f^ which corresponds to the expected 
weak ferromagnetic correlation in the canted spin liquid. 

One of the most striking features of the results is the presence of a gap in the spin 
excitations, without which the data could not be fit with an integer power law dispersion 
relation. Corroborating evidence for the existence of such a gap can be found in the earlier 
a.c. susceptibility data which showed that Xac decreases with decreasing temperature in the 
spin liquid phase [[T|,p3[. Further physical justification for the existence of the gap can be 
seen in the inset to figure [|, which shows that that the fitting parameter A(H) is linear 
in field as would be expected for an energy gap. We find that A(H) can be fit well by 
A(H) = A — /iexcH, and we associate \i exc = -0.16 \ib with a negative effective moment 
of the spin excitations which is much smaller than the moment of the S = 7/2 Gd spins. 
This finding is consistent with the fundamental magnetic unit of the spin liquid being a 
low-moment correlated cluster of spins in the spin liquid as has been suggested theoretically 
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The extrapolated zero-field gap, A /kB ~ 0.03K, could be attributable an intrinsic gap 
in the spin liquid or the anisotropy associated with the relatively large dipole interactions 
between the Gd spins. Our value of A /ks ~ 6 cw /70 is somewhat smaller than that 
predicted theoretically for the intrinsic gap in a kagome system [||, but it is of the same 
order. If we attribute Aq instead to a finite-size effect, we can estimate the correlation 



length of the excitations, Lq ~ ^4it 2 D/Aq ~ 40 A . This length scale is comparable to 
the ~ 100 A correlation length observed in recent zero-field neutron scattering studies ||. 
Interestingly, despite the first order transition between the spin liquid and the long-range- 
ordered antiferromagnetic phase at higher fields, the same functional form can fit C(T) in 
the long-range-ordered phase at low temperatures with a somewhat larger value of D, and 
A ~ 0. This suggests both that A originates from a finite-size effect in the spin liquid 
phase and that the short range antiferromagnetic order in the spin liquid phase is of the 
same sort as that in the long-range-ordered phase. 

In summary, we have investigated the low temperature excitations of a three dimensional 
spin liquid through its thermodynamic properties. We find that there is a crossover temper- 
ature below which the thermal fluctuations change character, which is apparently a generic 
feature of geometrically frustrated spin liquid states. We also find that we can describe the 
magnetization and heat capacity data by a simple model in which the excitations have a 
gapped quadratic dispersion relation analogous to ferromagnetic spin waves. These data 
provide the first thermodynamic evidence for the nature of the excitation spectrum in a 
three dimensional spin liquid phase pi ], and they indicate the need for new single-crystal 
inelastic neutron scattering investigations of these systems. 
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FIG. 1. The field dependence of static magnetization and the magnetic susceptibility dM/dH 
at 0.1 K. The inset shows the low-temperature phase diagram of GGG including the spin glass, 
spin liquid, and antiferromagnetic long-range ordered phases [3,9]. 
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FIG. 2. Specific heat (C) of GGG versus applied field at several temperatures, note 8C p /dT 
changes sign at ~ 0.235 K). 



8 



FIG. 3. The temperature dependence of C p and C p /T in the spin liquid regime in different 
applied fields. Note that all the curves cross at T 

— Tcross ~ 0.235 K. The inset shows the energy 
gap (A) obtained from fits to the low temperature C(T) data 
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FIG. 4. The temperature dependence of M (normalized to 0.06 K) in the spin liquid regime. 
The inset shows dM/dT versus temperature, notice the minimum at ~ 0.3 K corresponding to 
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FIG. 5. Left: Low-temperature M versus temperature at 0.4 and 0.6 T. The solid lines are fits 
to the data as described in the text. Right: Low-temperature C p versus temperature at 0.4 and 
0.5 T respectively. The solid lines are fits to the data as described in the text. 
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